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l ^ Abstract 

Skeletons of branching processes are denned as trees of lineages characterized by 
an appropriate signature of future reproduction success. In the supercritical 

Ql case a natural choice is to look for the lineages that survive forever [6]. In 

the critical case it was earlier suggested [7] to distinguish the particles with 
the total number of descendants exceeding a certain threshold. These two 
definitions lead to asymptotic representations of the skeletons as either pure 

i i birth process (in the slightly supercritical case) or critical birth-death processes 

__ (in the critical case conditioned on the total number of particles exceeding a 

high threshold value). The limit skeletons reveal typical survival scenarios for 

^O the underlying branching processes. 

o 

{-**«) In this paper we consider near-critical Bienayme-Galton- Watson processes 

^H- and define their skeletons using marking of particles. If marking is rare, 

such skeletons are approximated by birth and death processes which can 

be subcritical, critical or supercritical. We obtain the limit skeleton for a 

sequential mutation model J8] and compute the density distribution function 

p\. for the time to escape from extinction. 
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1. Introduction 

Imagine a population of viruses trying to establish itself in a new environment. 
Suppose the currently dominating type is nearly critical, in that its mean offspring 
number is close to one. One can think of two main factors which may lead to survival 
of this population: reproductive success or an advantageous mutation (a mutation 
producing new type of particles forming a strictly supercritical process). While a 
reproductive success is possible in the slightly supercritical case, 'survival due to an 
advantageous mutation' is the only way to escape extinction for a slightly subcritical 
branching system. 

The typical survival scenarios of such branching processes can be studied in terms 
of the so-called skeleton trees formed by lineages characterized by an appropriate 
signature of future reproduction success. In the supercritical case a natural choice is to 
look for the lineages that survive forever [6] . In the critical case it was earlier suggested 
[7] to distinguish the particles with the total number of descendants exceeding a certain 
threshold. These two definitions lead to asymptotic representations of the skeletons 
as either pure birth process (in the slightly supercritical case) or critical birth-death 
processes (in the critical case conditioned on the total number of particles exceeding a 
high threshold value). 

In this paper we suggest an alternative approach of defining a skeleton that relies 
on a random marking of the lineages in the family tree of a Bienaymc-Galton- Watson 
(BGW) process. The skeleton is then defined as the subtree formed by the infinite 
lineages together with the marked lineages. In Section [2] we describe the skeleton of 
infinite lineages and recall the result from [B] concerning a sequence of single type 
slightly-supercritical BGW process. It says that, conditioned on the event that the 
skeleton is not empty, the skeleton is approximated by the standard Yule process (a 
linear pure birth process). 

In Section|3]we consider an exactly critical BWG process with marking: each particle 
in the family tree is marked, independently of the others, with a small probability. 
Here the skeleton is formed by the lineages leading to a marked particle. We show 
that, conditioned on the event that it is not empty, the skeleton is approximated by a 
critical linear birth-death process. 
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In Section HI we state the main result of the paper concerning a sequence of near- 
critical BGW process with marked particles. The definition of skeleton is adapted in 
order to include both infinite lineages and marked lineages. The marking is also done in 
a more general way than in Section [3j Our main result states that, conditioned on the 
event that it is not empty, the skeleton is approximated by a linear birth-death process 
which can be either supercritical, critical or subcritical, depending on the parameters 



of the model. The proof of this main result, Theorem 4.1 is given in Section 8 



In Section [5] a decomposable two-type BGW process with irreversible mutations, 
starting from a single wild type individual, is studied. Each daughter of a wild type 
individual becomes a mutant, independently of the others, with a small probability. 
We look at this Binomial mutation model as a particular case of the processes treated 
in Section [4] by considering the following marking procedure: a wild type individual is 
marked if gives birth to at least one mutant daughter. 

Section \& deals with a sequential mutation model, considered in [5], for a viral 
population that escapes extinction due to a sequence of irreversible mutations that 
lead to a target type. It is assumed that mutations appear according to a Binomial 
mutation model and results from Section [5] are used to obtain the asymptotic shape of 
the limit skeleton. Finally, in Section 7 we address the important question of the time 
to escape from extinction in a sequential mutation model. By 'the time to escape from 
extinction' we mean the first generation where a particle of the target type appears. 
Due to the shape of the skeleton, we are able to show that the time to escape from 
extinction is asymptotically equal to the time of the first death occurring in the limit 
skeleton. An explicit formula for the density distribution function of the time to the 
first death in the limit skeleton is derived. 



2. Infinite lineages 

Consider a sequence of branching processes {Z m (n)}^L , m = 1,2, . . . with offspring 
distributions (p m (0),p m (l), . . .) and starting from one particle Z m (0) = 1. Assume that 



4 S. SAGITOV AND M.C. SERRA 

the processes are nearly critical with 

oo 

7^fcp m (fc) = l + e m , e m -> 0, m -> oo, (2.1) 

fc=i 

oo 

Y^ k(k — l)p m (k) — > a 2 , m — > oo, for some a > 0, (2.2) 

fc=2 

oo 

sup^ k 2 p m (k) ->• 0, n -> oo. (2.3) 



Condition (2.3) requires uniform integrability for the sequence of squared offspring 



numbers and implies that the following equality 

oo 

J2 s k p m (k) = 1 - (1 + e m )(l - s) + {a 2 /2 - R rn {s)){l - s) 2 (2.4) 

fc=0 

holds with R m (s) — $ uniformly in m as si" 1. 

A natural way of defining a skeleton for branching processes was proposed in [5], 
where such processes were considered with e m > for all m (slightly supercritical case) . 
Its survival probability, Q m , according to Lemma 3.3 in [6] satisfies the following well 
known approximation formula 

Q m ~ 26 m er~ 2 , m — > oo. (2-5) 

It is a well known fact that a supercritical branching process can be viewed as 
a two- type branching process, by distinguishing among particles with infinite line of 
descent and particles having finite number of descendants. If we concentrate only in 
the number of particles with infinite line of descent we arrive at the so-called skeleton 
process. Conditioning the supercritical process on non-extinction and focussing on 
infinite lineages we get a new sequence of supercritical branching process {X m (n)} n >o 



with X m (0) = 1 which under conditions (2.1 1, (2.2), and (2.3| is weakly approximated 



{X m (t/e m )} t > -> {*i(«)}t>o, m -> oo (2.6) 

by the Yule process, see Theorem 3.2 in [5] . Recall that the Yule process is a continuous 
time Markov branching process with particles living exponential times with mean 1 and, 



Skeletons 0} near- critical Bienayme-Galton-Watson branching processes 5 

at the moment of death, the particles are replaced by two new particles. The key part 



of the proof of (2.6) is to show, using (2.4), that 



P(X m (l) - l|X m (0) = 1) = 1 - e m + o(e m ), 
P(X m (l) = 2|X m (0) = 1) = e m + o(e m ). 

Then it remains to check the convergence of the generator of this Markov chain to the 
generator of the Yule process after the time is scaled accordingly. 

The limitation of this definition of a skeleton is that it has no direct extension to 
the critical or subcritical branching processes. Theorem 2.1 from [6 shows that in 
the critical case if the branching process is conditioned " on very late extinction" then 
the limiting skeleton (without any scaling) is a trivial discrete time process Y(n) = 1, 
n = 0,1,.... 

In this paper we suggest an alternative approach of defining a skeleton relying on a 
random marking of the lineages in the family tree. We start by studying in the next 
section a simple case of exactly critical reproduction. 

3. Critical branching processes with independently marked particles 

Consider a single type BGW process Z(n) such that its offspring distribution (p07.Pi! ■ • 
has mean X^fcLi kpk — 1 an d finite variance a 1 = *}^ =2 k(k — l)pk- Consider the 
corresponding family tree and suppose that each vertex in the tree is independently 
marked with a small probability fi m — > 0, as m — > 00. Any path connecting the root 
with a marked vertex will be considered as a part of the skeleton. Thus the skeleton 
is the subtree of the family tree formed by the skeleton paths. Adapting Proposition 
2.1 from [5] to this case one can show that for a given /i m the sequence {X m (n)} n >o of 
numbers of branches in the skeleton forms a BGW process. Next we find the conditional 
asymptotic structure of the skeleton. 

Let again Q m = ¥(X m (0) = 1) stand for the probability that the skeleton is not 
empty (at least one particle is marked). Due to the branching property we have 

1 - Qm = (1 " Mm)0(l - Qm), (3.1) 
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where (f>(s) = X}fcLiPfc s • Indeed, (3.1| simply says that the skeleton is empty if and 
only if the root is not marked and all the daughter subtrees, if any, have empty skeleton. 
Using the Taylor expansion of <fi around point 1 we get 



[I - Qm) = 1 - Qm + QW/2 + 0{Q 



2 ) 



and deduce from (3.1| 



Qm-v^V^- (3-2) 

Denoting by £ m the indicator of the event that the ancestral particle is marked we get 

E(r^s XmW ;X m (0) = l) =E(r^)E(/'»( 1 VE(r u s x "' (1, ;I m (0)=0) 
= (r Mm + 1 - Hm)<t>(sQ m + 1 - Q m ) - P(A m (0) = 0) 

implying that the offspring distribution of the skeleton particles satisfies 

v ( U x m (i)iy fri -s A -, l-0(l-Q,„(l-s)) u m 

E Ir^s mV >\X m (0) = 11=1 — (1 -r)— h o( v //^j 

1-r 



s + -y<9m(l - s) 2 - cry/2 /i m ■ — h o(^u^) 

(1 - cr\/2/i m )s + (T V / 2fi m f -r + -s 2 J + o{- s /Jx rn ). 



It follows that using r m = a~/2\x m we can write 



P(£ m = 0, A m (l) = l|X m (0) = 1) = 1 - T m + o{r m ), 
P(U = l,X m (l) = 0|X m (0) = 1) - r m /2, 
P(£ m = 0, X m (l) = 2|X m (0) = 1) - r m /2. 



(3.3) 



Comparing ( 3.3 ) to ( 2.7 ) we conclude that, if the original branching process produces 



at least one marked particle, there holds a weak convergence in the Skorokhod sense 



{X m (t/T m )} t > -> {Y Q . 5 (t)} t > 0l m 



oo. 



(3.4) 



Here the limit process is a continuous time Markov branching process with the critical 
binary splitting: 
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• particles live exponential times with parameter 1, 

• at the moment of its death each particle with probability 0.5 leaves no children 
and with probability 0.5 produces two children. 

The limit process, Y 5 (.), being a critical branching process will eventually go extinct. 



Relation (3.3 1 gives an enhanced interpretation of the limit skeleton (3.4). All 
marked particles appearing in the branching process can be associated with the tips 
of the limit skeleton. In particular the total number of the marked particles, W m , 
conditioned on W m > is asymptotically distributed as the total number of leaves, 
W , in the family tree of the limit skeleton. Due to the branching property we have 

W = 1 { „ =0} + (W + W") ■ 1 {V=2} , 

where v is the number of offspring of the initial particle in the skeleton, W' and W" 
are i.i.d. with W . In terms of the generating functions we get an equation E (s w ) — 
(s + [E (s w )} 2 )/2 leading to E (s w ) = 1 - JT- ~s. 

4. Main result 

In this section we combine and further develop the two approaches presented in 
Chapters [2] and [3] for a more general model. Consider a nearly critical reproduction 
law {p m ( k ) } ^L satisfying (2.1), (2.2), (2.3) and allowing for negative e m . Furthermore , 



assume that a particle with k offspring is marked with probability A m (k). The marking 
event may depend not only on the number of offspring but also on the whole daughter 
branching process. For example, the marking rule could be to mark all particles whose 
total number of descendent exceeds m [?]• Observe that for the marked near-critical 
BGW process the total probability for a particle to be marked is given by 



f'n 



22p m (k)A m (k). 



k=0 



Clearly, the case of Section [2] corresponds to the zero marking probability, /i m = 0, 
and the case of Section u| corresponds to e m = and A m (k) = [i m . 

Reconciling the two different definitions of a skeleton given in Sections [2] and [3] we 
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next introduce a new definition. 

Definition 4.1. For a given family tree of a BGW process with marking, the subtree 
formed by the marked lineages together with infinite lineages will be called the skeleton. 

Clearly, if /i m = the skeleton is formed only by the infinite lineages. If n m > 
any infinite lineage becomes marked and we can think that the skeleton is formed only 
by the marked lineages. 

In this paper we study the asymptotic behavior of the skeleton assuming 

H m -> 0, m -> oo, (4.1) 

restricting ourselves to the cases when the mean offspring number for the marked 
particles 



M m = ft,,?- ^2 k Pm(k)A m (k) 



fc=i 

satisfies 

limsupM m < co. (4.2) 

m—^oo 

By this we exclude such extreme situations as, for example, when A rn {k) is of order 
fi m k 2 for large k and ^k?p m (k) —> oo. Observe also that given (2.3| and (4.1 1 



J2Hk-l)p m {k)A m (k)^0, (4.3) 

which is obtained by using the inequality 

oo oo 

y]k(k- l)p m (k)A m (k) < n 2 [i m + 2J k 2 p m (k). 



k=2 



Let, as before, Q m stand for the probability that the skeleton is not empty. Now we 
can state our main result claiming that, conditioned on the event that the skeleton is 
not empty, a weak convergence of the following form holds 

{X m (t/T m )} t >o -» {Yx(t)} t >o, m4», (4.4) 
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for a convenient sequence (r m ) m >o and convenient A € [0,1], generalizing both (2.6 1 



and (3.4). Here for a given A G [0, 1] the limit process is a continuous time Markov 
branching process with binary splitting: 

• particles live exponential times with parameter 1, 

• at the moment of its death each particle with probability 1 — A leaves no children 
and with probability A produces two children. 

Remark. Importantly as with (3.4) in Section |3j by the claiming (4.4) we implicate 
that asymptotically there is one-to-one correspondence among the marked particles 
appearing in the branching process and the tips of the limit skeleton Y\(.). In particu- 



lar, in (2.6) the limit skeleton has no tips implying that under the corresponding time 



scale we can not expect observing marked particles in the branching process. 



Theorem 4.1. Under conditions (2.1), (2.2), (2.3), (4.1), (4.2) assuming that there 
exists a finite or infinite limit 



c = lim e m /y/JI^ 



(4.5) 



(i) ifc=oo, then (2.5) and (2.6) hold, 
(ii) if c € (—00,00), then 



c+ \fc? + 2er 2 



(iii) ifc= —00, then 



l^m / I ^m I , 



also (A A) holds with r m = |e m | and A = 0. 



(4.6) 



also (4.4) holds with r m = w 'u m v / c 2 + 2<r 2 and A = h + \ , „ c „ ^ 



According to Theorem |4.1| there are five different asymptotic regimes for the skeleton 
of a near-critical BGW process depending on how the deviation from the critical 
reproduction, e m , relates to the square root of the marking probability, ^/i ra : 

• in the supercritical case c = 00 with a negligible marking probability the limit 
skeleton is the Yule process which never dies out, 
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• in the supercritical case c G (0, oo) with a balanced marking probability the 
limit skeleton is a supercritical Markov branching process which dies out with 



probability , c 2 + „ a 2 — and survives forever with probability 



2 c 



v / c 2 +2er 2 +c l J v / c 2 +2er 2 +c' 

• if the reproduction law is very close to the purely critical, c — 0, then the limit 
skeleton is a critical Markov branching process which dies out with probability 
one although rather slowly, 

• in the subcritical case c G (— oo, 0) with a balanced marking probability the limit 
skeleton is a subcritical Markov branching process which dies out with probability 
one, 

• in the subcritical case c = — oo with a very small marking probability the limit 
skeleton is given by a single lineage that dies out after an exponential time. 



5. Binomial mutation model 



Here we present an important example of a marked branching process based on a 
decomposable two-type Galton- Watson process modeling a population of individuals 
with irreversible mutations. The two- type branching process stems from a single wild 
type individual which produces k offspring with probability q m (k). Suppose that each 
daughter of a wild type individual becomes a mutant with probability 7r m independently 
of other daughters. 

To introduce a marked BGW process we focus only on the wild type individuals 
and mark those wild type individuals who have at least one mutant daughter. The 
reproduction law for the marked branching process is given by the distribution for the 
number of wild type offspring: 

p m (k) = ^ q m (k + ( fc I ') (1 - T m ) V ro , (5.1) 

and the conditional marking probabilities A m {k) are computed using the following 
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relations obtained by splitting (5.1| in two parts 



p m {k){l - A m {k)) = q m {k)(l - 7r m ) fe 
P 



(k)A m (k) = J2 ^(k + ( * 7 ' ) C 1 " ^H, 



i=i 



k + l 



(5.2) 
(5.3) 



To ensure that one can use the results from previous section, we need conditions 



(2.1)-(2.3) to hold. Therefore we assume that the reproduction law with mutant 



offspring satisfies 

oo 

},kq m (k) = 1 + 7/mi Vm -> 0, 
fc=l 

oo 

^k(k- l)q m {k) -> cr 2 , 

fe=2 

oo 

sup N^ k 2 q m (k) — > 0, n — » oo, 

m , 

A;— n 

for some a £ (0, oo). Assume also 

7Tm — > 0, m — >• oo. 



(5.4) 
(5.5) 
(5.6) 



(5.7) 



Lemma 1. Conditions (|5_4|), ( |5.5| ), ( |5.6[ ), ( |5.7| ) zmj% ( |2.1| ), ( |2.2[ ) ; | |2.3| ), ( |4.l[ ), ( |4.2 | 



Mm ~ 7Tm, TO ^ OO, 

M m — > cr 2 , to — > oo, . 



(5.8) 
(5.9) 



Remark. Relation (5.9 1 has an interesting implication for our two-type branching 



process: the mean number of wild type siblings in a family with at least one mutant 
asymptotically equals the variance of the total offspring number. 



Proof. First observe that (5.3| entails a useful expression for the marking probability 

oo oo 

Mm = ^P m {k)A m {k) = Y, 9m(*)(l - (1 - Km?)- (5.10) 



fc=0 



fc=l 
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Clearly, (5.10) and (5.4) yield 



< Mm - 7T m (l + Tj m ) < TT m E k 2 q m (k), 



(5.11) 



fc=i 



and (5.8) follows from (5.4), (5.6) and (5.11) 



Next, due to (5.2) we have 



< J2 k ( k - ^mik) ~ ]T k(k - l)p m (k)(l - A m {k)) 

k=2 k=2 

CX) 

= Y / k (k-l)q m (k)(l-(l-TT, n ) k ) 

fc=2 

oo 

< n ir m + sup 2_. k Qj(k) 



for any n > 2. Letting here first m — > oo and then n — > oo, due to (5.5) and (5.6), we 
arrive at 



E fc(fe - l)p ro (fc)(l - A m (k)) -> a 2 . 



fc=2 



This together with (4.3) implies (2.2). Now, according to (5.3) we have 



oo oo 



li m M m ^m = E fc ^~] q m (k + Q I ) (1 - 7r m ) fc 7r, 



fe=i z=i 



fc + ; 



\*s * — i 



E«™Ci)ECi-0(;)(i-* 



v'-M-i 



'm; " m 



3=2 1 = 1 

CX) 

(1 - TTm)- X J2 3U - l)Pm(j)(l " A„(j)) 



J=2 

oo j 



+ ^ e ?m0*) Ec? - o ( , ) (i - ^^ 

j=2 Z=2 



-Z_Z-2 



From here we easily obtain (5.9), and therefore (4.2), using (4.3) and 



E *»o*) EC? - o k) a - ^r'i 2 < E^'C? ~ ^^c?). 



j=2 i=l 



i=2 
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In view of 



^ kp m (k)(l - A m (k)) = 1 + e m - jU m M„ 



fc=l 



we derive from (5.2) and (5.4) 



< Vm ~ £ m + fti.M m < 7T m ^ fc2 9m(fc)- 



fc=l 



Combining (4.2) and (5.12| we get (2.1) 



To prove (2.3) we turn to (5.1 1 and see that 



oo oo 



J2 k 2 p m (k) = J2 ^ Yl qm( - k + Z ) 



k=n 2=0 
oo j—n 



k + l 



k„l 



(1 - 7T m ) fe 7T 






< ^2j 2 Qm(j)- 

j=n 



Thus (2.3) is an immediate consequence of (5.6) 



(5.12) 



Corollary 5.1. For the binomial mutation model satisfying (5.4), (5.5), (5.6), (5.7) 



and r/ m /^/ir m — > c the statements (i), (ii), (iii) of Theorem 4-1 are valid after (e m , /i m ) 
are replaced by (j] m ,ir m ). 



6. The sequential mutation model 



Our next illustration of Theorem 4.1 deals with the sequential mutation model [5] 
for a viral population with irreversible mutations which escapes extinction as soon as 
a target type of viruses is produced. To simplify the discussion we focus mainly on the 
two-step mutation model, extending the one-step model from Section [5j 

Suppose we have a population of viruses stemming from a single virus which is able 
to reproduce and mutate giving rise to what we call intermediate type of viruses. The 
viruses of intermediate type reproduce according to a common law and by mutation 
generate a new type of viruses which we call the target type. The marking rule for 
the intermediate type is straightforward: we mark mothers with at least one daughter 
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of the target type. The wild type marking rule is a bit more complicated: we mark 
a mother which has at least one successful mutant daughter (that is a mutant, of the 
intermediate type, which has at least one marked descendant in the whole line of 
descent). 

We will assume that the reproduction laws and marking probabilities for both wild 



type and intermediate type branching processes satisfy conditions of type (5.4), (5.5), 
(5.6), (5.7) and are described by triplets (r] m , c 2 , w m ) and (^mjcfi'^m) respectively. 



Furthermore, we assume for simplicity that 

Km/Km — > a, m — > oo, for some a G [0, oo). 
According to Corollary|5.1[ if 



then 



fjm ~ C2y/TT m , m — > oo, for some c 2 G (—00, oo), 
a , c 2 + a/c 2 + 2aa\ 



u 2 

where Q m is the counterpart of Q m in Corollary 



5.1 



for the skeleton of the intermediate 



type. Clearly, the intermediate type is supercritical iff c 2 > and a > 0. Notice that 
with a = the skeleton of the intermediate type is the Yule process. The time scale 
intermediate type skeleton is given by 



T m = V^m V C i + 2acrf . (6.1) 

Our main interest in the two-step mutation model is of course the limit skeleton 
leading to the target type. Therefore, we want to apply Corollary |5.1| once again to 
the branching system with the probability of a successful mutation for the wild type 
viruses given by 

=. _„, A ^3/2 C 2 + y/c 2 , + 2C 



Lcta 



r2 



Assuming 

Tjm ~ c i n m> m ~^ °°j f° r some Cl € (— oo, oo) 
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we get 

implying 

cH<7 2 

The limit skeleton for the wild type is supercritical iff ci > 0. 
The time scale for the wild type type skeletons is given by 



s + Vc 2 + 2<7 2 


c 




ClCT 2 


^-2 


V C 2 




<7 


+ Vc 2 + 2aa\ 


Ci(T2 + y cferf 


+ 2cr 


2 (c 2 + 




^4 + 2aa\) 




IT 5 


0-2 





/c 2 + 2a 2 a 2 - 2 ( C2 + Y , c 2 + 2aa 2 ), 
which in the considered case is much slower than the time scale of the intermediate 



type (6.1). Thus the overall skeleton is given by the wild type skeleton, and the first 
death in the limit skeleton corresponds to the time of escape from extinction when the 
first virus of the target type appears. 

The above considered case is one of the many possible combination of reproduction- 
mutation regimes for the two-step mutation model. Without analyzing each of the 
remaining cases we just point out that there is a situation when both parts of the 
skeleton live on the same time scale. This is the case when 

Tt m ~ cot^, to — > oo, for some 7 G (1, 2), 
fj m ~ — fiir^ 1 , to — > 00, for some /3 <E (0, 00), 
T)m ~ Ci7r m , to — > 00, for some c\ € (—00,00). 

Here both time scales are of order l/7r m . 

Turning to the the sequential model with 6—1 intermediate steps before the target 



type, we extrapolate the formula (6.2) to Q m ~ const • 7T^~ . This prediction should 
be compared with the strictly subcritical case where one expects Q m ~ const • 7r^ , see 
Theorem 7.1 in [8]. 
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7. Time to escape 

In the framework of the sequential mutation model of Section [6] it is crucial to be 
able to describe the time to escape from extinction. 

Let T m be the time until the first marked particle is observed and put 



Q(t) = lim p(T m > 



t 



v / /i m (c 2 + 2(7 2 ) 



According to Theorem 4.1 given that the limit (4.5) is finite, Q{t) = P(T > t) is 
the tail probability of the time T to the first death in the limit skeleton Y\(-) with 
A = | + | , 2 C + 2 :. The branching property ofY\(-) says that 

r = X + min(T , ,T")-l {v=2}) (7.1) 

where L is the exponential life length with mean one, v is the number of offspring of 



the initial particle, T" and T" are i.i.d with T. Due to the branching property (7.1 1 



Q(t) = P(L >t) + P(T >t,L<t) 



-t i \ / r&i 



e't + X Q z {t~u)e- U du. 
Jo 



It follows 



t 



e t Q{t) = l + \ / Q\u)e u du. 
Jo 

Differentiation over t yields a simple differential equation 

Q'{t) + Q(t) = \Q 2 (t), Q(0) = 1 

giving Q{t) = 1/(A + (1 - A)e*). Thus 

t \ / c / c 



7^/ V Vc 2 + 2a 2 v Vc 2 + 2a 2 



'- 



t v / c 2 +2o- 2 



Skeletons of near- critical Bienayme-Galton-Watson branching processes 
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and we conclude that the scaled time T m J[i rn has the limit density distribution function 

1{y/<? + 2ct 2 - c ) e *Vc*+2 CT * 



m = 



1 



v / c 2 +2<t 2 

In particular, with c = we get 



Vc 2 +2ct 



= ") e *Vc 2 +2(T : 



r, t>0. 



</>(*) = 2V2a 2 e ta ^(l + e t °^)~ 



(7.2) 




Figure 1: Three density curves given by formula ( 7.2 \ with a — 1: subcritical case c = —0,5 



(solid line), critical case c = (dashed line), supercritical case c = 0.5 (dotted line). 



We illustrate the asymptotic density function, ip, by Figure [T] In the supercritical 
case the density curve reaches its a maximum value at 



*E 



1 



In 1 



2c 



Vc 2 + 2a 2 '" V ' Vc 2 + 2a 2 - c 
making the most likely value for the time to escape T m to be around 



Tm ~ V /Uc 2 + 2a 2 ) ln V 1 + 7^T2^ - c 
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8. Proof of Theorem UTTl 

Our proof relies on the properties of the probability generating function 

oo 

f m (r,s) = Y,Pm(k)s k (l-A m (k)+A m (k)r) 

k=0 

jointly characterizing the marking status of a particle (through r) and its offspring 
number (through s). 



Lemma 8.1. Given conditions (2.1|, (2.2|, (2.3), and any sequence u m € (0,1) such 



that u m — > 0, as m — > oo, the following decomposition holds 

1 - fm(r, 1 - u m ) = n m (l - r) + u m (l + e m - /i m M m (l - r)) - u m a^ m /2, 
with 07 m ~ * (j2 as m ^ co uniformly over r <G [0, 1]. 



Proof. This follows from a Taylor expansion around point (r, 1) 

f m (r, 1 - Um) = f m (r, 1) - u m ^(r, 1) + ^ ( ^%(r, 1) + i? m (0 



9s v ' ; 2 V 5s 
where 

oo 
fm(r, 1) = 5^Pm(*)(l - A,„ (A) + A m (fc)r) = 1 - Mm(l - f), 

fe=0 

df °° 

-|=.(r, 1) = J2 k Pm(k){l - A m (k) + A m (k)r) = 1 + e m - \x m M m {\ - r), 

fc=i 
-^(r, 1) = J2 fc ( fc - !W( fc )(! - A m{k) + A m (k)r), 

k=2 

and 

CO 

R m (r) = £fc(fc - l) Pm (k)(l - A m (k) + A m (k)r)(l - O^ 2 ), 

k=2 



Skeletons of near- critical Bienayme-Galton- Watson branching processes 
for some 9 G (1 — u m , 1). Indeed, since for any n > 3 

oo 

Rm{r) <Y, k2 Pm(k){l - (1 - 
k=2 

oo 

< u m n 3 + 2J k 2 p m (k), 
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\fc-2\ 



fc— n 



condition (2.3 1 implies i? m (r) — > 0. It remains to apply (2.2) and (4.3) 



The skeleton is empty if the initial particle is not marked and all her children produce 
empty skeletons 

Using Lemma 8.1, with r = and u m = Q m , we obtain a quadratic equation 



^Q.m^Vm ^^CmX^m l^m^m ) ^^m U, 



entailing 



Q,. 



£m - MmMn + \/(fm - H, n M m ) 2 + 2a^ m fi r , 



where <7o, m — > c Once again applying Lemma 8.1 now to the right hand side of 
E[r^s x ™W;X m (0) = 1] = E[r^s x ^] - P[X m (0) = 0] 

= f m (r, sQ m + 1 - Q m ) - (1 - Q m ) 

with u m = 1 - Q m (l - s), we get 

E[r«-s*-«|X m (0) = 1] = (1 + e m - Q m a 2 ) s + Q m s 2 a 2 /2 + ^„Q m V (8.2) 

- MmQm 1 - Cm + Qm^ 1 j2 + o(Q m ) + 0(/X m ). 



Now we are ready to verify the statements of Theorem |4.1| case by case 



Case (i) If c = oo, then is \i m = o(e m ) and (8.1) yields (2.5). Furthermore, 
A^rnQm 1 = o{e m ) and (pT2J) gives 



E[r ? ms x m (i)| Xm(0) = i] = (i _ em ) s + €r 



(8.3) 
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as desired. 



Case (ii) If e m ~ Cy/fi m for c e (—00, +00), then (8.1 1 yields (4.6) as we can neglect 



the terms involving /j, m M m . Furthermore, /i m Q r 



/n» 



c+^fc+c 



and the last terms in 



(8.2) are negligible 



■MmQm 1 _ £ ™ + QmCT 2 /2 + o(Q m ) + 0{^ m ) 

y/JI^cr 2 , , , c + Vc 2 + 2cr 2 



: + yjc + CO 2 



CyfjTr, 



+ o(Va^) =o(v / M^)- 



We conclude 



E[r .? ms x m( i)| Xm(0) = !] = ^ _ ^^ c 2 + 2^2) s 



ffl^y/c 2 + 2er 2 



1 1 



2 2 V c 2 + 2a 2 



1 1 



+ o( v / /^T) 



,2 2Vc 2 + 2a 2 , 

as desired. 

Case (iii) If c = -co, then a m :— e m — fi m M m takes negative values for sufficiently 



large m and we derive from (8.1) 



^ (J 0,ml J "m 



a a,m yy a m + ^ a 0,m^ 



m, ^m 



I En 



as stated. It remains to see that (8.2) gives 



E[r«-s*">«|X m (0) = 1] = (1 - |e m |)s + |e m |r + o(e m ), 



since ^ m Qm ~ \ e m\, Mm = o(|e m |) and Q m = o(|e m |). 



Acknowledgements 

The research of Serik Sagitov was supported by the Bank of Sweden Tercente- 
nary Foundation and the Swedish Research Council grant 621-2010-5623. The re- 
search of Maria Conceigao Serra was financed by FEDER Funds through "Programa 
Operacional Factores de Compctitividade - COMPETE" and by Portuguese Funds 
through FCT - "Fundacao para a Ciencia e a Tecnologia", within the Project Est- 



Skeletons of near- critical Bienayme-Galton-Watson branching processes 21 

C/MAT/UI0013/2011. 

References 

[1] Athreya, K. and Ney, P. (1972). Branching processes, John Wiley & Sons, London-New York- 
Sydney. 

[2] Ethier, S.N. AND Kurtz, T.G. (1986) Markov Processes: Characterisation and Convergence, 
John Wiley & Sons, London-New York-Sydney. 

[3] HACCOU P., Jagers P., Vatutin V.A. (2005). Branching Processes: Variation, Growth and 
Extinction of Populations, Cambridge University Press, Cambridge. 

[4] IwASA, Y., MlCHOR, F., NOWAK, M.A. (2003). Evolutionary dynamics of escape from biomedical 
intervention. Proc. Roy. Soc. London B 270, 2573-2578. 

[5] IWASA, Y., MlCHOR, F., NOWAK, M.A. (2004). Evolutionary dynamics of invasion and escape. 
Journal of Theoretical Biology 226, 205-214. 

[6] O'Connell, N. (1993) Yule process approximation of the skeleton of a branching process. J. 
Appl. Prob. 41, 725-729 

[7] Sagitov, S. (1997) Limit skeleton for critical Crump-Modc-Jagcrs branching processes. In 
Classical and modern branching processes, IMA Vol. Math. Appl., 84, pp. 295-303, Springer, 
New York. 

[8] SAGITOV, S. and Serra, M.C. (2009) Multitype Bienayme-Galton-Watson processes escaping 
exctinction. Adv. Appl. Prob. 41, 225-246 

[9] Sevastyanov, B. A. (1971) Vetvyashchiesya protsessy (in Russian), Nauka, Moscow. 
Also available in German: 
Sewastjanow B. A. (1974) Verzweigungsprozesse, Akademie-Verlag, Berlin. 

[10] Taib, Z. (1992) Branching Processes and Neutral Evolution, Lecture notes in Biomathcmatics 
93, Springer, Berlin. 



